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In this article we present a closed analytical description for few-cycle, focused electromagnetic 
pulses of arbitrary duration and carrier-envelope-phase (CEP). Because of the vectorial character of 
light, not all thinkable one-dimensional (ID) shapes for the transverse electric field or vector potential 
can be realized as finite energy three-dimensional (3D) structures. We cope with this problem by 
using a second potential, which is defined as a primitive to the vector potential. This allows to 
construct fully consistent 3D wave-packet solutions for the Maxwell equations, given a solution of 
the scalar wave equation. The wave equation is solved for ultrashort, Gaussian and related pulses 
in paraxial approximation. The solution is given in a closed and numerically convenient form, based 
on the complex error function. All results undergo thorough numerical testing, validating their 
correctness and accuracy. A reliable and accurate representation of few-cycle pulses is e.g. crucial 
for analytical and numerical theory of vacuum particle acceleration. 

I. INTRODUCTION 

Recent developments in laser technology [1] resulted in ultrashort electromagnetic pulses, which may contain only 
a few optical cycles and can be focused down to a single wavelength leading to the so-called A 3 -regime. In addition, 
there was a tremendous experimental progress in carrier-envelope-phase (CEP) control of these ultrashort laser pulses 
allowing to synthesize almost arbitrary pulse shapes [3, H, 0] . Applications for these well controlled laser pulses range 
from coherent attosecond control B to high-gradient electron acceleration [f| 0, H, El and generation of ultrashort 
coherent X-ray flashes [Io|, EH El, US]. At the same time, numerical studies demonstrate the importance of correct 
analytical description of laser pulses in vacuum 0,01. It was shown, that even weakly inaccurate solutions of Maxwell 
equations can lead to largely erroneous results when applied blindly to, e.g., direct particle acceleration by the laser 
fields. Clearly, there is a demand for accurate analytical description of these pulses. 

As in the case of strongly focused pulses the vectorial character of light becomes crucial for few-cycle pulses. 
There is a significant interdependence between the pulse shape and the polarization that requires careful analysis. We 
demonstrate that not all field structures conceivable in ID models can be realized as finite energy, localized 3D wave 
packets. Commonly used approximations are consistent for a certain choice of the carrier envelope phase (CEP) only. 
When trying to construct laser pulses starting from a given shape for the transverse field component, one easily ends 
up with a pulse inconsistent in the 3D geometry. Our method to construct consistent 3D electromagnetic structures 
is valid for arbitrary CEP cases. 

Following the work of Porras [H, EH], our approach to the wave equation uses the analytic signal p0|. We also 
consider the particularly interesting case of a radially polarized laser pulse [H, EH- Because of its strong and 
purely longitudinal field component on- axis, the radially polarized pulse may become an important tool for electron 
acceleration 0, H, We provide the proper analytical solutions in a simple manner, which is particular convenient 
for use in numerical simulations. 

Finally, we let the solutions undergo some accurate numerical tests. Any significant errors in the solution would 
show up while they are propagated by the field solver. Compared to more conventional approximations, the new pulse 
description decreases electromagnetic artefacts drastically at the pulse initialization stage, and the self-consistent 
development of the pulse fields agrees with the analytical description in cases, where more conventional approximations 
fail. 



II. SECOND POTENTIAL REPRESENTATION 



An electromagnetic pulse can be represented by its four-potential A a = ((/>, A), where each component satisfies the 
vacuum wave equation DA a = 0. We use the Lorenz gauge d a A a =0 and further set the scalar potential to zero = 0, 
which can be done in vacuum. Then, the fields are written as E = —c~ 1 d t A and B = V x A. We are interested in 
finite energy pulse-like structures, so that A is required to be a localized function: \A\ — > for r — > oo. It is easy to 
see that the pulse potential is uniquely defined now, since each change of A generates measurable electric or magnetic 
fields. Because of <j) = 0, the Lorenz gauge coincides with the Coulomb gauge in vacuum V • A = 0. 

In laser physics, it is common to choose an analytical solution to the wave equation for the main, i.e. transverse 
components of the pulse. Then the longitudinal component can be determined: A z = Vj_ • A±dz. However, this 
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Figure 1: Logarithm of the square transverse magnetic field log 10 {[B y / Bo) 2 ) in the y = plane resulting from an ultrashort, 
linearly polarized, (a) sine- and (b) cosine-phased Gaussian potential. The pulse duration is cr = 0.5 A and focal spot width 
a = 2 A. 

integral can yield a non- vanishing longitudinal potential component far from the actual pulse region, since commonly 
used solutions to the wave equation do not satisfy the condition 

/oo 
V_l • A ± dz = (1) 
-oo 

for ultrashort pulses. Note that this component is not meaningless but will cause non-zero longitudinal electric and 
transverse magnetic fields, an example of which is shown in Fig. [J}). These fields have a small amplitude of the order 
O ((c/cjct) 2 ) , but since they extend infinitely along the beam axis, they contain an infinite amount of energy. 

To get a realistic finite energy pulse, our choice of the transverse vector potential is restricted by Eq. (pQ). This 
is a fundamental difference to the ID case, where such a restriction on the wave form does not exist. Before the 
consequences of this restriction are discussed in detail, let us introduce the second potential \l>, which enables us to 
describe a reasonable set of realistic pulse structures in a more convenient way: 

A = Vxf (2) 

Of course, each component of \l> has to satisfy the wave equation A^i — -^d^i = 0. Then, the wave equation for A 
and the Coulomb gauge readily follow from Eq. ([2]) and there are no restrictions like (pQ) on the choice of the second 
potential components. 

For a large class of laser pulses, it is convenient to choose e z \l> = 0. Then Eq. ([2]) becomes A z = —X7± • (e z x \l>) 
and = e 2 x 9 2 f . In the near-monochromatic case d z ~ — iu/c we get e z x \I> « icA±/uu, hence in the long 
pulse limit the transverse second potential (TSP) e z x \l> is, except for a constant factor, identical with the transverse 
components of the vector potential. 

Let us point out that each electromagnetic pulse in vacuum can be represented by the TSP. This is easily seen, 
since the second potential can be obtained by taking the integral e z x \l> = J z °° A^dz' from an arbitrary vector 
potential in Coulomb gauge. To make the representation unique, we require the condition |\l>| — > for |r| — > oo in 
the half-space x > 0, so that is unambiguously given by the just mentioned integral. Then, for a vast class of 
laser pulses including all linearly, circularly and radially polarized modes, \l/ will vanish at infinity in all directions. 
To understand this, we write the integral condition (pQ) in terms of the TSP: 



V ± -(e z x <Z\ z= _ oo )=0 



(3) 
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In general, it possesses non-trivial solutions, corresponding to structures, where the fields produced by \l/ at z — > — oo 
vanish, but \I> itself does not. However, we look at important special cases. For linear polarization (^ x = 0), Eq. ([3]) 
has obviously none but the trivial solution. Thus, all finite energy, linearly polarized pulses can be represented by 
a localized TSP. The same is true for circular polarization, which we define by Sk y = i^ x , and the components 
are assumed to be analytical functions. Eq. ([3]) then only has solutions of the type ^(x + iy), so there is no non- 
trivial solution that fulfills the boundary condition. Another interesting structure is the radially polarized pulse 
(e z x \l> = /(rj_, z,t)r±), and again there is no non-trivial finite energy solution to (0 obeying its symmetry. 

Summarizing, for each realistic linearly, circularly or radially polarized pulse, the condition e z x ^\ z= _ 00 = 

A±dz = on the transverse vector potential must be fulfilled. Commonly used analytical wave packet solutions 
strictly fulfill this condition only for a certain choice of the CEP. Take, for instance, a pulse with a Gaussian longitudinal 
profile, A y (x = 0, y = 0, z, t = 0) = aoexp [— (z/cr) 2 ] cos(27rz/A), then \^ y (z = — oo)| /(aoA) = 7r _1//2 r exp (— r 2 ) ^ 
0, meaning that some field components would extend to infinity as in Fig.Q}). Unlike the cosine-phased Gaussian, the 
sine-phased potential may in principle be assigned to A y , as the integral vanishes here (see also Fig. Hk). Assigning 
an exact Gaussian profile to the transverse electric field instead does not help. Indeed, by calculating the second 
potential of such a pulse it can be proven that this structure does not exist for any phase. 

Other than for linearly, circularly or radially polarized pulses, for the azimuthally polarized pulse e z x \l> = 
f(r±,z,t)e ip the TSP representation not have to vanish as z — > — oo. Notice that this pulse does not produce 
any longitudinal field at all. Anyway, it can neatly be represented by the longitudinal component e z • \l> instead of 
the transverse components of the second potential. 



III. SCALAR WAVE EQUATION 



Now we come to the solution of the scalar wave equation on \l>. The Fourier transform in time and the transverse 
directions yields d 2 ^ = — ((co/c) 2 — k±)S&, with the solution 



§(x,k ± ,c^) = § 2 _ Q exp ^-**y (^) - k lJ ( 4 ) 

Next, the focal spot profile ^{z = 0) has to be chosen in a physically reasonable way. Besides the very common linearly 
polarized Gaussian mode TEMoo, we also treat the very interesting radially polarised Hermite-Gaussian mode TMoi. 
Further we will tackle the 2D solutions, which differ in some factors from the 3D ones. Knowing the scalar solution 
for the TEMoo mode, we can easily construct circularly polarized Gaussian pulses by setting ^ y = i^ x , and moreover 
azimuthally polarized pulses by assigning the same term to the longitudinal component e z • \l/ of the second potential. 

Let ip(t) be the pulse time dependence at the center of the focal spot and its Fourier transform. The focal 

spot size may depend on the frequency: 



(e z x \l> J = exp - ^ 



r, TM 01 
TEMoo 



Now, the transverse direction Fourier transform of ([5]) is inserted into the solution of the wave equation ([4]). Since 

~ 2 
_L 



we consider pulses propagating mainly in one direction {(?k\ <C a; 2 ), we expand the square root in a Taylor series 



and neglect fourth order terms: y (oo/c) 2 — k 2 ^ ~ uo/c — ck\/(2u), performing the paraxial approximation. If now 

<j ex 1/y/oJ is assumed, so that the Rayleigh length zr\ = cocr 2 /2c is constant for all frequencies, we are able to carry 
out the inverse Fourier transform analytically and obtain the solution 

V ' ^ V Q J V c 2cq) \e x TEM 00 V } 

Here g = 1 for a linear (g = 0.5 in 2D) and g — 2 for a RP (g = 1.5 in 2D) laser and q — z + izr\ is the confocal 
parameter, ip is a complex representation of the time dependence of the pulse. Note, that since t' = t — z/c — r 2 L /2cg, 
it is generally a complex number. Choosing naively ip(t') = exp (— t /2 /r 2 + i(cjot f + </>)) yields a solution diverging 
for big r± as 0(exp(r^)). Instead the analytic signal [Tt} should be used, as suggested by Porras [15]. The analytic 
signal is the complex representation of a real signal without negative frequency components. The analytic signal 
representation of the Gaussian pulse g(t') = exp(— t /2 /r 2 ) cos(uuot + <p) is calculated to be 
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Figure 2: Direct comparison of a two dimensional cut through the complete short-pulse solution ([6]) using (a) the "naive" choice 
ip(t f ) = exp (— t /2 /r 2 + i(u;ot' + 0)) and (b) the analytic signal Eq. J7J). Pulse parameters are: cr = 0.5 A, a — 2 A, ct = —10 A 
(before focus) 



wherein r = and w (z) = exp(— z 2 )erfc(— iz) is the complex error (or Faddeeva) function (20|. The analytic signal 

of the Gaussian pulse was also calculated in [15], but Eq. (27) from [l5| disagrees with (0 for Im(^) 7^ 0. Analytical 
and numerical tests show, that ([7j) is the correct solution. To illustrate the meaning of ([7j) and its difference to the 
naive choice ip(t') — exp [—t ,2 /r 2 + i(uj$t' + (/>)), consider Fig. [2] In the region near the optical axis, Im(t') is small 
and the solutions agree quite well. However, for bigger r± the naive solution (b) diverges, while the analytic signal 
shows a proper beam-like behaviour and vanishes. 

IV. NUMERICAL TESTING 

The equations ([2]), ([6]) and (0 together form the key to accurate analytical and numerical representation of ultra- 
short few- and even single-cycle electromagnetic pulses. One important application of them is the use in numerical 
simulations and we will conclude this paper by showing their superiority to more conventional representations for this 
application, thereby checking the correctness and the accuracy of the analytical results obtained so far. 

We use the particle in cell (PIC) code VLPL [21] . To begin with, the fields are initialised inside the VLPL simulation 
grid. Then they are propagated using a standard algorithm on the Yee-mesh. Finally, it is verified if the numerically 
propagated pulse still agrees with the analytical term, overall or in some key parameters, and furthermore, if unphysical 
static fields remain at the place, where the pulse was initialized. 

First we check the correctness of Eqs. (0 and 0, proving their superiority to a conventional representation often 
used for numerical simulations. The "conventional", or separable form is a simple product of a monochromatic, 
transversely Gaussian beam with a Gaussian temporal profile. Circularly polarized pulses are used, so that the shape 
can well be seen in the intensity plots, with a duration of cr = A and a focal spot width of a = 2 A. The pulses are 
focused over a distance of 50A inside the simulation. In Fig. [3] the initial condition and the numerically propagated 
solution at the focal spot are shown. While the product approximation shows strong asymmetric deformation in the 
focus and the focal field does not reach its specified value of clq = 1, the proper short pulse solution is nearly perfectly 
symmetric in the focus and reaches the desired maximum. The presented solution Eqs. (|6]) and ([7j) is clearly superior 
to the simple product approach. 
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Figure 3: Focusing properties of different approximations of the wave equation, evaluated in a 2D version of the PIC code 
VLPL. (a) and (b) show the intensity distribution as it is initialised inside the code and (c) and (d) show the propagated 
solution at the focal spot, (a) and (c) use the CW paraxial solution multiplied with a temporal profile, (b) and (d) the correct 
short pulse solution. The circularly polarized laser pulses are Gaussian both in space and time with duration cr — A and width 
a 2A. 



Now we come to the second potential representation. The alternative to its use is to assign arbitrary wave equation 
solutions to the transverse components of the vector potential and then make some kind of approximation for the 
longitudinal part. The simplest possibility is to fully neglect the longitudinal field component, but for strongly focused 
pulses a somewhat more reasonable approximation can be reached by choosing A z = (c/iuoq)V ±_ • A_l, what we will 
call the quasi-monochromatic approximation, because it relies on d z ~ —iuoq/c. Using one of these approximations, 
the initial pulse structure has finite energy, and because of the energy conserving property of the field propagator 
algorithm, the pulse will be forced to self-organize into a consistent structure of finite extension. While this happens, 
"virtual charges" are left behind, unphysical static fields, which may make the concerned regions in the simulation 
domain unusable for further computations. We want to see, if this problem can be cured by the use of the second 
potential representation. 

When employing the second potential representation for an ultrashort pulse, the additional derivative will slightly 
alter the pulse shape. As shown before, this is inevitable, since a Gaussian shaped linearly polarized vector potential 
can exist as an independent structure in vacuum only if it is sine-phased. To make the TSP represented pulse 
comparable to the conventional one, it is necessary to take care for the frequency shift caused by the z-derivative, 
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Figure 4: Logarithm of the transverse ((a) and (b)) and longitudinal ((c) and (d)) square electric field log 10 [{Ei/ Eq) 2 ) , i £ 
{x, z) after a propagation time ct = 10 A. (a) and (c) depict the conventional, (b) and (d) the TSP Gaussian. Pulse parameters 
are cr — 0.3 A and zri — 4ty 2 A. For the second potential, a corrected frequency uj = ujo — y/2/r was used to take account for 
the frequency shift caused by the additional derivative. 

which can be estimated as cc; e fr = c^o + V%/t. 

Fig. [H shows the pulses after a short propagation distance in the PIC simulation box. Firstly one observes that, 
despite of the very short duration, the moving pulse structures (right half of the images) appear very similar in 
the conventional and the TSP version. Secondly one notices, that the conventional pulse leaves behind a significant 
amount of virtual charge fields in the initialisation region (left half of the images), having both a longitudinal and a 
transverse component. This undesired phenomenon can neatly be suppressed by the use of the second potential, seen 
in the figure. 

The last test we want to present in this paper concerns the longitudinal field component on the optical axis of a 
radially polarized pulse, which is of particular interest for vacuum electron acceleration @, !, ©|. The pulse length 
used was cr = 0.5A, and the focal spot size a = 2A. Again, the conventional approach corresponds to a near- 
monochromatic approximation A x — zV_l • A±/co, so as to generate a finite pulse structure from the given transverse 
field components. The phase (j) in Eq. ([7j) was chosen as <j) = for the second potential representation and as (j) = 0.5 7r 
for the conventional representation, so that the pulses are actually comparable. Initially, the longitudinal field nearly 
agrees for both representations, since the first term of Eq. ([7j), which is the dominating one, is the same in both cases. 

In figure [5] we present the longitudinal field of the pulse after it has left its "virtual charges" behind. The fields 
initialized using the conventional representation already differ significantly from the analytical description, whereas 
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Figure 5: The longitudinal pulse field in PIC at the time c(t — to) = 4 A after initialization, compared to the analytic solution. 
Once the TSP was used for initialization (blue line) and once the near-monochromatic approximation for A z . The differences 
in the analytical representations are too small to be seen in the diagram, so that both are represented by one curve (dashed 
black line). Pulse parameters are: cr = 0.5A, a = 2 A, 



the TSP represented fields agree almost perfectly. This will be crucial e.g. when PIC simulations are to be compared 
with other analytical or semi-numerical calculations, where the field is inserted analytically and is not self-consistently 
propagated. Without an exact and reliable pulse representation, such a comparison is hardly possible. 



V. CONCLUSION 



In this paper we have given a comprehensive guide to the mathematical representation of ultrashort, Gaussian 
and related, electromagnetic pulses. The wave equation for an ultrashort Gaussian pulse has been solved in paraxial 
approximation. Further it has been shown, that the vectorial character of light has a stringent influence on its 
field structure for ultrashort pulses. To pay regard to this, ultrashort pulses should be represented by their second 
potential, using the transverse components for linear, circular or radial polarization and the longitudinal component 
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for azimuthal polarization. Numerical tests proove that the solutions flawlessly work in the regime of ultrashort, 
moderately strong (a > A) focused pulses. 
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